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ABSTRACT 

We study the distribution of the durations of gamma-ray bursts (GRBs) in the unified model of short and 
long GRBs recently proposed by Yamazaki, loka, and Nakamura. Monte Carlo simulations show clear bimodal 
distributions, with lognormal-like shapes for both short and long GRBs, in a power-law as well as a Gaussian 
angular distribution of the subjets. We find that the bimodality comes from the existence of the discrete emission 
regions (subjets or patchy shells) in the GRB jet. To explain other temporal properties of short and long GRBs, 
the subjet parameters should depend on the angle in the whole jet. 
Subject headings: gamma rays: bursts — gamma rays: theory 



1. EVITRODUCTION 

The durations of gamma-ray bursts (GRBs) observed by 
BATSE show a bimodal distribution, which has led to a classi- 
fication of GRBs into two groups: bursts with Tgo durations < 
2 s are called short GRBs, and those with durations > 2 s are 
called long GRBs (Kouveliotou et al. 1993; McBreen et aT] 
Il994h . If Tqo directly reflects the active time of the progen- 
itor of the GRB, different origins of short and long bursts are 
implied, such that the former arise from the binary neutron 
star mergers while the latter arise from the collapse of mas- 
sive stars (e.g. Meszaros 2002; Zhang & Meszaros 2004). 

The short and long bursts roughly consist of 25% and 
75%, respectively, of the total BATSE GRB popula- 
tion. We should regard these fractions as comparable, 
considering possible instrumental effects on the statis- 
tics. If these two phenomena arise from essentially 
different origins, the similar number of events is just 
by chance. However, some observations have suggested 
that the short GRBs a re similar to the long GRBs (e.g., 
[ Germany et al.' '200 0(: E^azzati. Ramire z-Ruiz. & Ghisellini' 

boon iN akar & Pi ran) |2Q0 2|: iLamb et al. 2003; 

iGhirlanda. G hisellini. & Celotti 2004). Motivated by 
these facts, .Yamazaki. loka. & Nakamura ( ,2004b) proposed 
a unified model of short and long GRBs, even including 
X-ray flashes (XRFs) and X-ray-rich GRBs, and showed 
that it is possible to attribute the apparent differences of 
the light curves and spectra of these four kinds of events to 
the different viewing angles of the same GRB jet. This is a 
counter-argument against the current standard scenario of the 
origins of short and long GRBs. 

In this paper, we perform Monte Carlo simulations to show 
that our unified model naturally leads the bimodal distribu- 
tion of the Tgo durations of GRBs. The paper is organized as 
follows. In §12] we begin with a brief review of our unified 
model of short and long GRBs. The Tgo duration distribution 
is calculated in §|3l Sectionl^is devoted to discussions. 

2. UNIFIED MODEL OF SHORT AND LONG GRBS 

We briefly describe our unified model of short and long 
GRBs (for details, see Yamaz aki. loka. & Nakamura 2004t|). 
We assume that the GRB jet is not uniform but made up of 
multiple subjets, and that each subjet causes a spike in the 



observed light curv e. This is an extreme case of an inhomo - 
geneous jet model ('Nakamura' "200^; 'Kumar & Piran" '20001). 
Let us consider a subjet with the opening half-angle A6'sub 
moving with Lorentz factor 7, observed from the viewing 
angle 0y. Because of relativistic effects, the subjet emission 
becomes dim and soft w hen 9^ is larger than ~ A6'sub + 7~' 
(!loka & Nakamural l200ll) . The effective angular size of its 
emission region is 7r(A0sub + 7~')^, which is larger than the 
geometrical size of 7rA6'sub^- For the multiple subjet case, 
the crucial parameter is the multiplicity (n,) of the effective 
emission regions along a line of sight . If many subjets point 
toward us (i.e., ^ 1) the event looks like a long GRB, while 
if a single subjet points toward us (i.e., n., = 1) the event looks 
like a short GRB. 

Below we give a typical set of parameters for the tempo- 
ral and spatial configurations of the GRB jet to demonstrate 
which type of event is observed depending on n^. We suppose 
that A^tot subjets are launched from the central engine of the 
GRB randomly in time and directions and that the whole jet 
consists of these subjets. We introduce a spherical coordinate 
system {r^d,ip) in the central engine frame, where the origin 
is the location of the central engine, and = is the axis of 
the whole jet. The axis of the yth subjet (J = 1 , • • • , A^tot) is 
denoted by (?9'^' , ip'''), while the direction of the observer is de- 
noted by (i^obs, </5obs)- We suppose that the jth subjet departs at 
time fjgp from the central engine and emits at radius r = r''' and 

time f = t^^f.p + r"^ / [3"^c. The departure time of each subjet f^^p 
is assumed to be homogeneously random between f = and 
t = fdur, where fdur is the active time of the central engine mea- 
sured in its own frame and is set to fdur = 20 s. The emission 
model for each subjet is the same as the uniform jet model 
in Yamazak i. loka. & Nakamura (2003a). For simplicity, all 
the subjets are assumed to have the same intrinsic luminosity 
and opening half-angle A6'|,^(j = 0.02 rad, and the other prop- 



erties are 7" 



100, 



3 X 10^^ cm 



-1, P'i; 



-2.5, 



and ^hv Q = 500 keV for all j. The opening half-angle of 
the whole jet is set to A^tot = 0.3 rad. We randomly spread 
A^tot = 350 subjets following the angular distribution function 
of the subjets as 



He, < I? < 



(1) 
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where i?/, = A6'tot - A6'sub and dc 



= 0.02 rad (see also 
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Fig. 1 . — Angular distribution of Mot = 350 subjets confined in the whole 
GRB jet in our simulation. Each subjet is located according to the power-law 
distril)ution function of eq.Q- The whole jet has an opening half-angle of 
A9to( = 0.3 rad. The subjets have the same intrinsic luminosity and opening 
half-angles ASs^i, = 0.02 rad, and the other properties are 7 = 100, r = 3 X 
10'^ cm, «B = -1 , /3b = —2.5 and 'yhu'o = 500 keV. The effective angular 
size of the subjets are represented by the solid circles, while the whole jet is 
represented by the dashed circles. The examples of lines of sight A and B are 
shown in the figure, while C is located at (-0.04 rad, 0.04 rad) and D is close 
to the center of the whole jet. 



iRossi. Lazzati & Reesll200l IZhang & M eszaros"2002^. Fig- 
ure Q] shows an example of the angular distribution of the 
effective emission regions of the subjets in our calculation. 
Most of the subjets are concentrated near the 1!) = axis (i.e., 
the multiplicity in the center ~ 100). For our adopted pa- 
rameters, isolated subjets exist near the edge of the whole jet, 
and there are some directions in which no subjet is launched. 

Figure|2]shows examples of the observed light curves in the 
50-300 keV band, each of which corresponds to the lines of 
sight A, B, C, and D in Figure^ The coordinate (i^obs, '/'obs) 
of C is (-0.04 rad, 0.04 rad), and D is close to the center of 
the whole jet. If many subjets point in the direction of the line 
of sight, such as in the cases of C («., = 15) and D (n^ = 97), 
we see a spiky temporal structure. In the case of B (wj = 2), 
the event consists of the distinct emission episodes. These 
are identified as long GRBs. If only one subjet points to- 
ward us, like in the case of A (n^ =1), the contributions to 
the observed light curve from the other subjets are negligible 
because of relativistic beaming effect, so that the observed 
gamma-ray fluence and duration are both about a hundredth 
of the typical values of long GRBs. These are quite similar to 
the characteristics of short GRBs. In addition, when the line 
of sight is away from any effective subjet regions (i.e., = 0), 
the soft and dim prompt emission is observed because of rel- 
ativistic Doppler effect and beaming effe ct, which is iden- 
tified as an XRF or an X-ray-rich GRB ("loka & Nakamura' 
pgOl; Yamazaki. loka. & Nakamura 2002, 2003a, 2004a b: 
mimazaki. Yonetoku. & Nakamural2003bl) . 

3. DISTRIBUTION OF Tgo DURATION 

We perform Monte Carlo simulations to show that our uni- 
fied model can explain the observed bimodal distribution of 
Tgo durations of GRBs. We fix the subjets' configuration 
as in Figure ^ We vary only the line of sight of the ob- 
server and calculate the Tgo duration for each observer in 
the 50-300 keV band. We generate 2000 lines of sight with 
< i^obs < 0.35 rad according to the probability distribution 
of sin i?obs iii9obs dipobs- We then select only hard events, whose 
observ ed hardness ratio is 5( 2-30 keV)/5(30-400 keV) < 
10""^ dSakamoto et alJl2004l) . The other soft events are clas- 
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Fig. 2. — Observed light curves in the 50-300 keV band for the lines of 
sight shown in Fig.0 A with rts = 1 (upper left), B with n,, = 2 (upper right), 
C with ris = 15 (lower left), and D with «j = 97 (lower right). The sources are 
located at z = 1 . The Tgo durations are 0.25 s for A, 14. 1 s for B, 25.4 s for C, 
and 37.8 s for D. 
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Fig. 3. — Distribution of multiplicity for the angular distribution of the 
subjets of Fig. Q The dashed line represents the analytical estimate of the 
n,"^ line (see text). 



sified as XRFs or X-ray-rich GRBs, which are observed when 
all subjets are viewed off-axis. 

Figure |3] shows the distribution of «.s in our simu- 
lation. The multiplicity «s is roughly proportional to 
n(i?obs, V'ob.s)- Then the distribution of n., is given by P(ns) oc 
sm('do\,^){d'do\,s/dns) ^ (Fig.|3 dashed line). We first con- 
sider the Tgo distribution in the case in which the redshifts of 
all the sources are fixed at z = 1 for simplicity. The result 
is shown in Figure One can see a bimodal distribution of 
T90 clearly. Which type of burst is observed, long or short, 
depends on Wj, and the distribution of is unimodal. Then 
why does the distribution of the duration become bimodal? 
The reason for the scarcity of the events for 1 < Tgo < 10 s 
is as follows. Let us first consider the event with «s = 1. In 
this case the Tgo duration does not vary significantly around 
^ 0.25 s when 6',, < A6'sub, which is determined by the angular 
spreading time of a subjet. As the viewing angle increases, T90 
increases (loka & Nakamura 2001). When ^i, > A^sub + T ^ 
however, the emission becomes soft and dim, so tha t the event 
will not be d etected as a GR B ( Yamazaki. loka. & Nak amural 
2002, 2003at lYamazaki. Yo netoku. & Nakamura 2003bj). The 
T90 takes a maximum value of ^ 0.75 s when 6^. ^ A0sub+7~'- 
We confirm that = 1 for almost all Tgo < 1 s events. Next let 
us consider the = 2 case. The example of the light curve for 
this case is Figurel^?, and the Tgo is 14.1 s. The Tgo duration 
is roughly given by the interval between the arrival times of 
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Fig. 4. — TgQ duration distribution in the 50-300 keV band of hard events 
with observed fluence ratio 5(2-30 keV)/S(30-400 keV) < 10"" ^ The jet 
model is the power-law. All sources are located at z = 1 . The dashe d line 
represents the analytical formula for the long GRBs, given by eq. IA2I . 





100 




80 




60 


1 




3 




z 


40 




20 









0.01 



0.1 



1 

To 



10 



100 1000 



Fig. 5. — Same as FigRlbut the source redshifts are varied according to the 
cosmic star formation rate(see text for details). Both short and long GRBs 
look like lognormal distributions. 



two pulses. Since the two pulses arrive sometime in the range 
< Tibs < Tdur, where Tdui is the active time of the central en- 
gine measured in the observer's frame, Tdur = (1 +z)fdur = 40 s, 
the mean interval is 40/3=13.3 s. This means that the duration 
of the tis = 2 event is much longer than that for Wj = 1 . For 
n.s > 3, the mean duration is longer than 13.3 s. The typical 
example is Figure|2^ for «, = 15, with Tgo = 25.4 s. This is the 
reason we have few events for 1 < 7yo < 10 s. The maximum 
value of Tgo is ~ T^ur- For the long bursts, the distribution 
function of Tgo durations can be derived from a simple proba- 
bility argument (see the AppendixlAlfor details). The dashed 
line in Figure l^represents the analytical formula of equation 
iA2\ . On the other hand, the distribution function of the short 
bursts seems to be too complicated to calculate analytically, 
since it sensitively depends on the jet configurations, such as 
the angular distribution and the intrinsic properties of the sub- 
jets. 

The ratio of events of the short GRBs and the long 
GRBs is about 2 : 5, which can be explained as follows 
£famazaki, loka, & Nakamura 2004b). The event rate of the 
long GRBs is in proportion to the effective angular size of the 
central core t?e.e//' ^(0.15 rad)^, where n, > 2. The event rate 
of the short GRBs is in proportion to M( A^sub + 7~')^, where 
M is the number of isolated subjets in the envelope of the 
core and M ^ 10 in our present case. Then the ratio of event 
rates of the short and long GRBs becomes M(A6'sub + 7~')^ : 
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Fig. 6.- 
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events with observed fluence ratio 5(2-30 keV)/S(30 - 400 keV) < Ur"^. 
We calculate the ON time duration as the time during which the emission 
is larger than 10% of the peak flux. The subjet distribution is given by the 
power-law form. The source redshifts are varied according to the cosmic star 
formation rate. 



W-^2:5. 

In reality, we should take into account the source redshift 
distribution. We assume that the rate of GRBs is in proportion 
to the cosmic star formation rate. We adopt the model SF2 in 
|R)rciani & Madau (2001), in which we take the standard cos- 
mological parameters of Hm = 0.3 and Ha = 0.7. Figure |5] 
shows the result. The distribution is again clearly bimodal, 
and the shapes of the short and long GRBs look like lognor- 
mal distributions. The ratio of the number of short and long 
GRBs is about 2 : 5 in this case as well. The dispersion of the 
lognormal-like distribution seems relatively small compared 
to the observations. This is ascribed to simple modeling in 
this paper We fix the jet configuration and use the same in- 
trinsic properties of the subjets. If we vary fdur for each source 
and 7"' for each subjet randomly, for example, the dispersion 
of lognormal-like Tgo duration distribution will increase from 
the general argument that the dispersion of the lognormal dis- 
tribution increases with th e increase of the nu mber of the as- 
sociated random variables Jloka & Nakamural2002 ). In more 
realistic modeling the observed dispersion will be reproduced. 

4. DISCUSSION 

We have investigated the Tgo duration distribution of GRBs 
under the unified model of short and long GRBs proposed 
by Yamazaki, loka, & Nakamura (2004b), and found that the 
model can reproduce the bimodal distribution observed by 
BATSE. In our model, the crucial parameter is the multiplicity 
(ris) of the subjets in the direction of the observer The dura- 
tion of an «s = 1 burst is determined by the angular spreading 
time of one subjet emission, while that of an n, > 2 burst is de- 
termined by the time interval between the observed first pulse 
and the last one. These two different time scales naturally lead 
a division of the burst Jgo durations into the short and long 
ones. We also performed a similar calculation for a Gaussian 
distribution, n(-d,(p) = ricCxp [-(i9/?9c)^/2] , and found that the 
790 duration distribution is bimodal in the same way as for the 
power- law subjet model. 

Let us make another compari son of our model with BATSE 
data. iMitrofanov et alJ Jl998h have computed the distribu- 
tion of the observed pulse number (denoted by Hp in their 
paper) and found that it is unimodal. If the Hp distribution 
were compared with the distribution, our model might be 
compatible with the observations, although some long GRBs 



4 



Toma, Yamazaki, & Nakamura 



80 
60 

B 

20 

- 




0.01 0.1 1 10 100 

On-time duration (sec) 

Fig. 7. — Same as Fig.|S| but the intrinsic luminosity of each subjet is 
assumed to be A'^' =Ao(= constant) fori? < 0.15 radand A"' = Ao(i?/0.15)"'' 
for 0.15 rad < i9 < where Aq is in arbitrary unit. Then we only take the 
events with peak flux larger than 3 X 10""* of the maximum peak flux that has 
appeai'ed in the calculation. The dashed line represents n, > 2 events. 
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Fig. 8. — Same as Fig. |5| but the emission radius r''' of each sub- 
jet is assumed to be r'^' = 3 X lO'^ cm for i5 < 0.15 rad and r''' = 3 X 
10''*(i?/0.15)^ cm for 0.15 rad < i? < The source redshifts ai'e fixed 
as z = 1. The solid line represents all the events, while the dashed line repre- 
sents > 2 events. 



are identified as np = 1 events. They also derive the distri- 
bution of the ON time duration — defined as the time during 
which the emission is larger than 40% of the peak flux — 
and found it bimodal. Furthermore, they argue that the mean 
pulse widths of short and long GRBs are different. On the 
other hand, we computes the on time duration distribution in 
the context of our theoretical model and found it unimodal 
(see Fig. |6}, which is expected since the pulse widths are 
almost the same. However, there are several observational 
implications that the distances to short GRBs detect ed with 
BATSE are smaller than those of long GR Bs (e.g., iTavanil 
[T998':'Ghir landa. Ghisemni. & Celo tti'fooT). although this is 
controversial. Then the observed pulse widths for short and 
long GRBs might be different because of the redshift factor 
To give an example, let us assume that the intrinsic luminosity 
of each subjet in the core region of the whole jet is larger than 
that in the periphery of the whole jet and count only the GRB 



events with peak flux larger than 3 x 10 of the maximum 
peak flux in our simulation. The result is shown in Figure 
in which we find that the effect of the peak flux cut off con- 
tributes to the bimodality of the on time duration distribution. 

At present, the observationally inferred bimodality of on 
time duration is not explained in our current model, in which 
all the subjets have the same intrinsic luminosity, the same 
opening half-angle, the same gamma factor, the same emis- 
sion radius, and so on. This is an extreme modeling for sim- 
ple calculation. In reality, they may depend on the off-axis 
angle in the whole jet; so may the pulse widths. Furthermore, 
iNakar & Piran (2002) investigated the pulse widths of GRBs 
using 2 ms time resolution and report that short GRBs also 
consist of several pulses. This can be incorporated into our 
model by assuming that a subjet radiates successive emissions 
rather than one instantaneous emission. Then the pulse width 
with 64 ms resolution (which is used in Mitrofano v et alJ 
Il998h will be determined by the active time of the subjet. 
If the pulse widths from the subjets in the central part are 
larger than those in the periphery, the bimodality of on time 
duration distribution can be explained. For example, we as- 
sume that the emission radius r'" is larger for the core region 
than for the periphery. Figure |8l is the result, which shows 
the bimodal-like distribution. Therefore, as we show in two 
examples (Figs. and |8}, some modifications of our model 
contribute to the bimodality of the on time duration, so that 
the current observed on time duration distribution is not in- 
consistent with our model. We hope that in the future more 
sophisticated modeling will reproduce the observed on time 
duration distribution. 

It has commonly been said that the observed bimodal dis- 
tribution of the Tgo durations of BATSE bursts shows the dif- 
ferent origins of short and long GRBs. However, the bimodal 
distribution is also available as a natural consequence of our 
unified model of short and long GRBs. The clear prediction 
of our unified model is that short GRBs are associated with 
energetic supernovae (SNe), since the associati on of long du- 
ration GRB s with SNe is strongly suggested (iGala ma et al 
[1998; Stane k et alJ2003t [Hiorth et al. 2003; Delia Vall e et al 
l2003i) . Indeed, one of the short GRBs shows possible associ- 
ation with a SN (Germanv et al. 2000). Even if the SNe are 
not identified with short GRBs because of some observational 
reasons, we predict that the spatial distribution of short GRBs 
in host galaxies should be similar to that of the long GRBs. 
Another prediction is that short GRBs have the same total ki- 
netic energies as lo ng GRBs, which might be confirmed by 
radio calorimetry dBerger et alJ2003l) . 
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14204024 (T. N.). 



APPENDIX 

ANALYTICAL ESTIMATE OF THE INTRINSIC Tgo DISTRIBUTION OF THE LONG BURSTS 

In this Appendix we derive the analytical distribution function of the Tgo durations of the long GRBs when all sources are 
assumed to be at z = 1. At first we consider for a given «.,(> 2). Each subjet causes one pulse, whose shape is a J-function for 
simplicity. In the present case the arrival time of the pulse from each subjet is random in the range < rob,s < Tdur- For a given 



Bimodal distribution of GRBs 



Tgo, the first pulse is required to arrive within Tdui - Tgo. The arrival time of the last pulse is determined as the time Tgo after the 
first pulse. The rest of the pulses are required to arrive within the range of Tgo. Thus, the probability function of Tgo for a fixed 
is approximately given by 

PnXTgo)dTgQ = n,(n,-l) — [7^] (Al) 



dur \ dur / ^ dur 

2 



For the power-law angular distribution of the subjets the distribution function of Wj is proportional to , so that we get 

D/T -2d (79o/7dui)+[l-(79()/7dur)]log[l-(ryo/rdur)] 

P(T<)o)dT<)o cx } n, P„XT9o)dTt)o = — ■ (A2) 

^ ^9()/^dur ^90 

The distribution function of n, for the Gaussian angular distribution of the subjets can be obtained in a similar way. 
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